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SUMMARY 


A simplified physical model is constructed which simulates the viscous 
crossflow in a fluid layer near the slots at a fixed streamwise location 
in a slotted wind tunnel. For low to moderate Reynolds numbers, numerical 
solutions of the two-dimensional, incompressible Navier-Stokes equations 
in stream function and vorticity, which govern the model flow, are 
obtained. Fairly general slot geometry is incorporated by means of the 
Thompson-Thames-Mastin transformation. An approximate factorization 
scheme with cyclic acceleration parameters is employed to solve a 
finite difference analog of the stream function equation. The vor- 
ticity equation is numerically solved with a modified version of the 
classical alternating direction implicit (ADI) scheme. Although no 
quantitative assessment of solution accuracy can be made, numerical 
results for variations in incremental wall pressure around the slat 
are at least qualitatively similar to some experimental results of 
Berndt and Sorenson [21] . 
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INTRODUCTION 


Longitudinal slots are often employed in wind-tunnel test sections 
to reduce adverse wall effects. However, the capability to correct 
wind-tunnel data in order to account for boundary- induced flow distor- 
tions is still rather limited. Analyses resulting in formulations of 
the slotted wall boundary condition [1-6] are usually based upon 
inviscid representations of the flow. On the other hand, it is 
conjectured by some that viscous effects of an unknown and complicated 
nature are present; therefore, results from the inviscid models become 
suspect [7]. The inviscid theory has been rather well developed and 
some experimental verification of flow models has been attempted 
[1-11]. Nevertheless, alleged inconsistencies in trying to apply 
theoretical results in practical flows often result in placing 
reliance upon empirical methods in the design and use of slotted 


test sections. 


The purpose of the present investigation is to provide some insight 
into the effects of viscosity on slot flows. The approach employed is 
a viscous crossflow analysis in a transverse plane. Suppose it is 
assumed than an aerofoil extends the breadth of the tunnel. By reason- 
ing analogous to the slender body theory, it is argued classically 
[2-6, 11] that the inviscid perturbation crossflow at a fixed longi- 
tudinal station in a test section with multiple slots above and/or 
below the model can be approximated in a fluid layer near the slots 
by solution of the two-dimensional velocity potential equation. The 
present treatment is a viscous generalization of such flow models, 
through use of the two-dimensional, time-dependent, incompressible 
Navier-Stokes equations and numerical solution techniques. 

This approach does not satisfactorily model the flow, particularly 
since streamwise (down- tunnel) boundary- layer growth and interaction 
is not accounted for. However, some insight into Reynolds number 
effects and flow separation in the vicinity of the slot should be 
provided. Due to the complicated geometry involved, ordinary 
finite difference numerical methods are expected to encounter 
difficulty with treatment of boundary conditions unless coupled 
with transformation techniques. The Thompson-Thames-Mastin (TTM) 
transformation [12] is an aid in attaining more accurate results 
for general shapes. By this method the geometric restrictions 
which are inherent in classical flow models due to idealizations 
in slat shape, e.g., infinitely thin slats and slots with sidewalls 
(infinite depth), can be avoided. 
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NUMERICAL ANALYSIS 


The Flow Domain 

The computational domain and a cross-sectional view of the slot 
geometry for the two-dimensional problem are provided by figure 1. 

Flow symmetry with respect to centers of slot and slat is assumed. 

The problem is driven by slat-generated vorticity and by a prescribed 
inflow at the top of the computational domain, on the tunnel side of 
the flow. We remark that prescribing the inflow also prescribes mass 
flow rate at outflow, although the strict form of velocity distribu- 
tion at the outflow can be altered from what it is at inflow. Boundary 
conditions for the case of identically prescribed mass flux distribu- 
tions at inflow and outflow are provided by figure 1(a). 

Numerical Algorithm 

In order to simplify the numerical solution procedure, the domain 
of figure 1(a) is mapped onto a rectangular region by means of the TTM 
transformation. A typical finite difference mesh generated for the 
problem using the Thompson grid-generation package is given by figure 
1(b). Taking the transformed coordinates to be x and y, the two- 
dimensional, time-dependent, incompressible Navier-Stokes equations 
in dimensionless stream function- vorticity form for the transformed 
problem are now given: 
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Here the functions of x*y, denoted by a, 0, y, a, t* J, are transforma- 
tion coefficients; R is the Reynolds number; £ is vorticity; if) is 
stream function; and the quantities u,v are related to transformed 
velocity components. Exact expressions for u*v and the transformation 
coefficients are given in reference 13. 

The problem now is the numerical solution of equations (l)-and (2) 
on a rectangular region. Whatever algorithm is applied to equation (1), 
the overall efficiency is determined by the manner in which equation (2) 
is solved* which solution can occupy as much as 90 percent of the total 
time required [14]. Classically* this difficulty has been alleviated 
through use of fast Poisson solvers applied to equation (2) . Such 
methods are inapplicable in the case of variable transformation coeffi- 
cients. Therefore* in the present study the method chosen for solving 
equation (2) is an approximate factorization scheme [15] which amounts 
to the use of the alternating direction implicit technique (ADI) with 
cyclic, accelerated convergence parameters [16]. This method has 
proved superior to the popularly used successive line overrelaxation 
techniques for the equations of transonic potential flow [15] . 

Let 


C = 
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- 1 


( 3 ) 


be the correction at stage n to an estimate °f the stream function 
at a lattice point* and let be the residual when equation (2) is 
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approximated using second order, centered difference derivative 
approximations, and is substituted. The stream function itera- 
tion is given by 

(S - a<S )(S - y5 )C = SwR (4) 

Here $ xx >&yy are second order, centered difference operators and equa- 
tion (4) results in tridiagonal systems of linear equations to be 
solved in the customary alternating direction fashion. With oj = 2, 
the parameter S is cyclically varied to speed convergence by reducing 
error- frequency components in a band determined by the maximum and 
minimum magnitudes of the transformation coefficients a,y. At the 
interior nodes, the values of the quantities u,v and the transforma- 
tion coefficients a, 3, y, a, t, J are obtained by centered differ- 
encing of the expressions given for these variables in reference 
13. 

Once having obtained the velocity-related quantities u and v 
from a stream function solution of equation (2) , the vorticity can 
be updated by solving equation (1). For efficiency, the time step 
employed should be the maximal restricted only by requirements of 
accuracy, and the unconditional linear stability of the classical 
ADI schemes lends itself to satisfaction of this criterion. We 
have therefore chosen to solve equation (1) using Briley's imple- 
mentation [14] for the ADI scheme of Douglas [17]. 

The ADI scheme is applied by splitting the time step into two 
equal parts. Over the first half-step the following equation approxi- 
mates equation (1) : 
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The equation for the second half-step is 
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The approximations (5a, b) become second order accurate in space and time 

when u* . , v* . and the mixed derivatives Z are evaluated at time 
i,3 i,3 xy 

4+1/2' This i s accomplished by half-step predictions of the form 
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Second order accurate space differencing is applied to 


5 xy 
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For equations (5a, b), boundary values of vorticity are prescribed zero 
except along the wall (which is transformed to the right boundary of the 
computational plane) . Here a vorticity distribution at the new step is 
computed using the full step prediction equation (7) . 
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Upon solving equations (5a, b) and updating the stream function through 
equation (4), a corrected vorticity distribution is calculated. This 
process is repeated several times, to converge wall vorticity prior to 
proceeding to the next step. On the first correction the new vorticity 
distribution at the wall is obtained by using second order accurate back- 
ward differences in equation (2) wherever necessary. Thereafter the 
Israeli correction [18] 


r m+l m /a_ 
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( 8 ) 


is applied to iteratively induce the no-slip condition. Here the deriva- 
tive at the wall is evaluated third order accurate using backward 

o X 

differences, since near separation points convergence of the iteration was 
found to be sensitive to this parameter. The parameter K is constant for 
all wall points, but is increased in magnitude as tangential velocity 
decreases. The proper choice of the sign of K that will induce no-slip 
is discussed in reference 19. 


COMPUTATIONAL RESULTS 

In this section preliminary numerical solutions are presented for 
viscous slot flow at Reynolds numbers of 100 and 1000, with a slot thick- 
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ness ratio and slot openness ratio (slot length to slot spacing) of 0.1. 
Openness and thickness ratios are realistic; however the fully rounded 
slot shape of figure 1 is not known to have been used by wind tunnel 
designers. Customarily, slat shapes have been nearly rectangular or 
trapezoidal, with small radii of curvature (as low as 1/160 of the slat 
length) produced by rounding sharp corners. 

Figure 2 exhibits comparisons between viscous and inviscid slot flow 
stream function and vorticity profiles. Figures 3 and 4 show pressure 
distribution around the slat, relative to plenum-side stagnation pressure. 
The quantity plotted is 


AP(y) = P(y) 


P 

s 



- 8£ y ) dy 


(9) 


where y is the transformed coordinate which varies along the wall and P^ 
is the plenum pressure at the center of the slat. The y-plot scale has 
been normalized to vary from zero to one around the slat, with 0.5 the 
coordinate of slot center. Figure 5 shows the character of the boundary 
layer across the center of the slot for the Re = 100 normal component of 
velocity. Values of -v are plotted. The inviscid normal velocity distri- 
bution across the slot is also shown. 

These results were obtained using a 51 by 51 grid, which is uniform 
in the transformed plane, but highly irregular in the physical plane, 
with order of magnitude variations in step size [see fig. 1(b)]. We 
remark that a boundary-layer packing function has been applied so as to 
concentrate the majority of points near the wall and in the vicinity of 
the slot. The scale of figure 2 is approximately actual size, or 1 . 0 x 
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8.2. The calculations were stopped when the maximum change in wall 
vorticity became the order of lO" 4 , with a maximum wall velocity value 
of magnitude 10“ 4 . 

Approximately 700 time steps using At = 0.0001 were required to 
converge the Re = 1000 calculation, with starting flow field initialized 
using the Re = 100 solution. The code required 153 k g core storage 
locations, of which approximately 31 k g were used to store metric coeffi- 
cients and transformation-related quantities. Typical Cyber 175 average 
CPU time for the advancing one time step is 2.2 x 10~ 4 seconds per step 
per computational node per model equation. Of course, to advance each 
time step requires a number of iterations to converge wall vorticity in 
order to yield an acceptably small level of wall velocity, and each 
vorticity iteration requires a few iterations (an average of four after 
impulsive start effects damp out) to converge the stream function field. 

No real comparisons of efficiency between algorithms can be made unless 
the same problem is solved by another method; however, it appears that 
Rubin’s modification [20] to the strongly implicit method of Stone could 
perhaps produce more efficient convergence at the expense of some addi- 
tional core storage. 

Flow calculations at Re = 5000 became unstable. This is probably due 
to insufficient mesh refinement. It is expected that, as the Reynolds 
number increases, thinning of the boundary layer near the lip of the 
slot in the vicinity of minimum wall pressure can quickly lead to 
problems with mesh refinement. It is not known at which point turbulent 
transition brought on by increasing Reynolds number might produce physical 
instability. 
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SIMILARITY TO EXPERIMENT 


Consider the calculated slat pressure distributions of figures 3 
and 4. As the flow enters the slot, a distinct pressure drop resulting 
from geometric flow acceleration is observed. A slight recompression 
precedes flow separation. For the case Re = 100, a pressure plateau 
occurs, followed by a gradual pressure recovery on the plenum side of 
the slat. However, the inference implied by the pressure data for Re = 
1000 is that, as the Reynolds number increases, viscous effects due to 
the separation bubble produce yet another pressure drop prior to the 
final recovery stage. 

It is expected that this second drop becomes increasingly sharp at 
higher Reynolds numbers, continuing the trend observed from figures 3 and 
4. Some support for this inference may be gleaned from experimental data 
of Berndt and Sorenson [21], which is reproduced in figure 6. 

Although slat shapes, slat thickness and openness ratios, and 
Reynolds numbers are somewhat dissimilar, comparison of figures 4 and 5 
(with pressure axis positive-direction reversed for figure 5) shows 
qualitative agreement for pressure data measured close to the aerofoil 
(position x = -22.5). 

DISCUSSION AND CONCLUSIONS 

In any numerical study the accuracy of results is subject to ques- 
tion, particularly as in cases such as this where theoretical results 
are not available for comparison. Grid refinement has not been attempted, 
although our feeling drawn from previous numerical experiences is that the 
grid is a little coarse in the vicinity of the slot, and possibly near 
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plenum stagnation. However, the qualitative similarity of numerical and 
experimental wall pressure distributions is certainly encouraging. Further- 
more, the results of figures 2 and 3 appear qualitatively reasonable. The 
flow separates in a region of adverse pressure gradient; the pressure then 
temporarily flattens in the separation region as the separation point 
(body distance equal to 0.54) is passed. The separation point is past 
the lip of the slot on the plenum side, and the circulation bubble is 
clearly defined. The flow does not appear to reattach; this may be due to 
numerical error, perhaps caused by the wall vorticity iteration. We 
conclude that overall the numerical algorithm has performed well, 
although its efficiency is hindered by the wall vorticity calculation. 

The approximate factorization algorithm for the stream function calcu- 
lation appears remarkably efficient, requiring very few iterations to 
convergence once the calculation damps out the impulsive start effects. 

Finally, we emphasize that this report represents only a cursory 
study of the viscous slot flow problem. Further efforts appear necessary, 
with much that can be done using the present code in terms of investi- 
gating the effects on the flow of slot geometry and Reynolds number. 
Moreover, three-dimensional models which incorporate the effects of 
turbulence and streamwise boundary layers are necessary for meaningful 
insight into the slot flow mechanism. 
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INFLOW: f = -VqX, $ = 0 


SLOT 

SYMMETRY 
\jt= 0 

5 = o 



Figure 1(a). Computational domain and boundary conditions. 




Figure 1(b). Coordinate system (alternate lines parallel to wall omitted) 
6x1 magnified, horizontal x-direction. 
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Figure 2. Slot flow profiles: viscous vs. inviscid 
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